ANALYSIS OF A MATHEMATICAL MODEL OF SYNTROPHIC BACTERIA 

IN A CHEMOSTAT 

TEWFIK SARI, MILED EL-HAJJI, JEROME HARMAND 

Abstract. A mathematical model involving a syntrophic relationship between two populations 
f^^ , of bacteria in a continuous culture is proposed. A detailed qualitative analysis is carried out. The 

f^ ■ local and global stability analysis of the equilibria are performed. We demonstrate, under general 

assumptions of monotonicity, relevant from an applied point of view, the asymptotic stability 
of the positive equilibrium point which corresponds to the coexistence of the two bacteria. A 
syntrophic relationship in the anaerobic digestion process is proposed as a real candidate for this 
^— ,. ■ model. 

(N 

^. 

Q ', 1. Introduction 

r^ . A synthrophic relationship between two organisms refers to a situation where the species exhibit 

"ti ' mutuahstism but where, at the opposite of what happens in a purely symbiotic relationship, one 

of the species can grow without the other. Such a situation can be mathematically formalized as 
follows. Assume that a first species denoted Xi grows on a substrate Si forming an intermediate 
product 5*2. This intermediate product is required by a second species X2 to grow. The limiting 
substrate of the second bacteria being the product of the first bioreaction, the second bacteria 
cannot grow if the first one is not present. 

Such interactions are quite common in nature: it is why a number of models have already been 
^\ . proposed in the literature. Katsuyama et al. [9- , proposed a model involving two mutualistic species 

10 ' for describing pesticide degradation, while a more general case is considered by Kreikenbohm and 

r->^ . Bohl [10 . Since mutualism involves generally species interacting through intermediate products, 

(^ ' other studies consider mutualistic relationships in food webs. For instance, Bratbak and Thingstad 

[2], or more recently, Aota and Nakajima [1^ considered the mutualism between phytoplankton and 
bacteria through the carbon excretion by the phytoplankton. A model studied by Freedman et 
al. [5] was proposed to explain the observed coexistence of such species. However, in the previous 
(«j ■ studies the models are very specific. In particular, the mathematical analyses of the models are 

i^ , realized for specific growth rates that are explicitely given (in most cases as Monod functions). 

C^ ' To extend the study of mutualism to more general systems, we have recently considered more 

general assumptions notably with respect to the growth rate functions considered in the models in 
using qualitative hypotheses, cf. [5] . Furthermore, it was assumed that the species Xi may be 
inhibited by the product 5*2 that it produces itself while the species X2 was simply limited by ^'2 . An 
example of such interactions was given by the anaerobic digestion in which mutualistic relationships 
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allow certain classes of bacteria to coexist. A mutualistic relation has been also considered in [3]. 
See [6] for another model of coexistence in the chemostat. 

In the actual paper, following [7], we revisit the model proposed in [5] in considering two main 
changes which significantly further extend the range of practical situations covered by the model. 
First, we assume that there is some S2 in the influent. In other terms, the limiting substrate ^2 
on which the species X2 grows is not only produced by the species Xi but is also available even 
if the species Xi is not present. The second modification of the model is that the second species 
is supposed to be inhibited by an excess of ^i, the limiting substrate on which the first species 
grows. To illustrate the usefulness of such extensions of the original model by El Hajji et al. [S], the 
biological interpretation of these hypotheses within the context of the anaerobic process is given in 
the appendix. 

The paper is organized as follows. In Section [2l we propose a modified system of four differential 
equations from the original model in [5]. The positive equilibria are determined and their local 
and global stability properties are established. In the case when the system has a unique positive 
equilibrium, the global asymptotic stability results are demonstrated through the Dulac's criterion 
that rules out the possibility of the existence of periodic solutions for the reduced planar system, the 
Poincare-Bendixon Theorem and the Butler-McGehee Lemma. Hence, in this case, for every positive 
initial conditions, the solutions converge to the positive equilibrium point which corresponds to the 
coexistence of the two bacterial species as observed in real processes. Simulations are presented 
in Section 21 an example of a syntrophic relationship is given in Section |5] as a candidate for this 
model. 



2. Mathematical model 

Let Si, Xi, S2 and X2 denote, respectively, the concentrations of the substrate, the first bacteria, 
the intermediate product, and the second bacteria present in the reactor at time t. We neglect 
all species-specific death rates and take into account the dilution rate only. Hence our model is 
described by the following system of ordinary differential equations : 

Si = D{Sr-Si)-hfii{Si,S2)Xi, 
Xi =. fii{Si,S2)Xi~DXi, 

S2 = D{S^" - S2) - k2tl2iSl, S2)X2 + kifIiiSi,S2)Xi , 
X2 = m{Sl,S2)X2~DX2. 



(1) 



^\ 



satisfies 



Where S*™ > denotes the input concentration of substrate, S*™ > denotes the input concentra- 
tion of the intermediate product and I? > is the dilution rate. 
Assume that the functional response of each species /ii,/i2 : I 

Al: /xi, /X2 : R+ -^ R+, of class C^ , 

A2: /ii(0,52)-0, M2(5i,0)=0, V(5i,52)eK^, 



A3: -91(^1,^2) >0, -^(5i,52)<0, y{Si,S2)& 

obi ob2 

A4: 1^(51,^2) <0, ^iSi,S2)>0, W{Si,S2)e 

obi ob2 



^l^ 
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Hypothesis A2 expresses that no growth can take place for species Xi without the substrate ^i and 
that the intermediate product ^2 is obhgate for the growth of species X2. Hypothesis A3 means 
that the growth of species Xi increases with the substrate 6*1 and it is inhibited by the intermediate 
product 5*2 that it produces. Hypothesis A4 means that the growth of species X2 increases with 
intermediate product S2 produced by species Xi while it is inhibited by the substrate 5*1 . Note 
that there is a syntrophic relationship between the two species. 

We first scale system ^ using the following change of variables and notations : 



fci 



ki 



si = —Si, xi = kiXi, 82^82, X2 ^ k2X2, s'^" ^ —S'i 



ki 
The dimensionless equations thus obtained are : 

DisY" - si) - fl{si,S2)xi 



ks 



^2 ~ '-'2 



(2) 



±2 



Where the functions /i , /2 



= /i(si,S2)xi - Dxi , 

= -D(s*2" - S2) - /2(si,S2)a;2 + /l(si,S2)xi , 

= /2(si,S2)a;2 - Dx2 . 
i^ — > M+ arc defined by 



fi{si,S2) ^ Hi{— 81,82) and /2(si, S2) = /i2(-;— si,S2). 

Ki fci 

Hypotheses A1-A4 satisfied by the functions ni and fi2 translate in the following assumptions of 
the functions /i and f2'. 

HI: /i, /2 : M^ ^ K+, of class C^ , 

H2: /i(0,S2)=0, /2(si,0) = 0, V(si,S2)eM^, 

H3: ^(si,S2)>0, ^{si,S2)<0, V(si,S2)eM^, 
O81 O82 

H4: |^(si,S2)<0, |^(si,S2)>0, V(si,S2)eK^. 
O81 O82 

M^, the closed non- negative cone in M^, is positively invariant under the solution map of system 
(IH). More precisely 

Proposition 1. For every initial condition in R5j_, the solution of system (0j has positive compo- 
nents and is positively bounded and thus is defined for every positive t. The set 

il ^ U8i,Xi,S2,X2) eM.^ : 81+ Xi ^ s"\ 82+X2 ^ Xi + 8^^> 

is a positive invariant attractor of all solutions of system (0j. 

Proof. The invariance of R^ is guaranteed by the fact that : 

i. si = ^ si = D si" > 0, 
ii. 82 = 0^82= D 4" + /i(si,0) xi > 0, 
iii. Xi = ^ Xi = for i = 1,2. 
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Next we have to prove that the solution is bounded. Let zi = si +xi, then ii = —D{zi 
which one deduces : 



(3) 



Sl 



{t)+Xi{t)^ 



Thus si{t) and xi{t) are positively bounded. Let Z2 
which one deduces: 



+ (si(0) + a:i(0)- 

S2 + X2 



-Dt 



Xi, then Z2 



-D{Z2 



s™) from 



from 



(4) 



S2{t) + X2{t) - xi{t) = 4" + (S2(0) + a;2(0) - xi(0) - s™)e 



-Dt 



Thus S2{t) and X2(t) are positively bounded. Hence, the solution is defined for all positive t. From 
([3]) and (jlj we deduce that the set fi is an invariant set which is an attractor. D 

3. Restriction on the plane 

The solutions of system ([2]) are exponentially convergent towards the set J7 and we are interested 
in the asymptotic behavior of these solutions. It is enough to restrict the study of the asymptotic 
behaviour of system ([2|) to il. In fact, thanks to Thieme's results [H], the asymptotic behaviour of 
the solutions of the restriction of ([2]) on 17 will be informative for the complete system, see Section 
[5l In this section we study the following reduced system which is simply the projection on the plane 
(a;i,a;2), of the restriction of system ([2|) on il. 

C xi = [$i(xi,a;2) -£)]xi, 
(5) < 

[ X2 = [^2(X1,X2) - D]X2. 

where 



$i(xi,a;2) = /i 



Xi, S2 



Xi 



X2) , $2(2;i,a;2) = /2 (sl" - Xi, S™ + Xi 



X2 



Thus, for ^ the state-vector (xi,X2) belongs to the following subset of the plane, see Fig. [T]: 
S= {(a;i,X2) eM+^ :0<xi <s\'',0<X2 < xi + 4"} . 



+ s% 




Figure 1. The set S 



The point F'^ = (0,0) is an equilibrium of ([5]). Besides this equilibrium point the system can 
have the following three types of equilibrium points. 

• Boundary equilibria F^ = (xi, 0), where xi = Si is a solution, if it exists, of equation 
(6) $i(a;i,0)=-D, 
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• Boundary equilibria F^ — (0, X2), where X2 = X2 is a solution, if it exists, of equation 

(7) $2(0,X2)=i?, 

• Positive equilibria F* = {x\,X2), where xi = x\^ X2 = X2 is a solution, if it exists, of the 
system of equations 

^ ' 1 ^2{XUX2)^D. 

We use the following notations 

^i-/iK",4"), i?2 = /2(sr,4")- 

The mapping xi n> $i(xi,0) is decreasing, and the mapping xi M- $2(2;!, 0) is increasing. If 
F>i > D2, there exists a unique real number ^1 satisfying $1(^1, 0) = $2(^1,0), since 

$1(0,0) = 1)1 > £'2 = $2(0,0), and$i(s5",0)=0<$2(sr,0). 

We denote by D3 &]D2,Di[ the unique real number (see Figure[2l right) such that: 

*i(a,o) = $2(6,0) = i?3. 



-I'2(xi,0) 



y = $2(2:1,0) 




y = <t>i{xi,0) 



Figure 2. Existence and uniqueness of xi. On the left, the case Di < D2: 
$2(^1,0) > D for all D < Di. On the right, the case Di > D2: $2(^1,0) > Z? if 
and only if D < D3 . 



The mapping X2 ^ $i(0, X2) is increasing, and the mapping X2 ^ $2(0, X2) is decreasing. Hence, 
if Di < D2, there exists a unique real number ^ satisfying $i(0,6) = $2(0,^2), since 

$2(0,0) ^D2>Di = $1(0,0), and $2(0,4") - < $i(0,4"). 

We denote by D4 g]-Di, 152 [ the unique real number (see Figure 121 right) such that: 

$1(0,6) = $2(0,6) = A. 
The nature of the trivial equilibrium point F^ is given in the following lemma. 

Lemma 1. If D > nia.x{Di,D2) then F^ is a stable node. // min(£'i, D2) < D < max(I?i,_D2) 
then F'^ is a saddle point. If D < imii(Di, D2) then F^ is an unstable node. 
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*l(0,i2) 




2/ = $1(0, 2:2) 



y = $2(0, a;2) 
X2sf 2:2 



Figure 3. Existence and uniqueness of X2- On the left, the case D2 < Di. 
$1(0, ^2) > D for aU D < Z?2. On the right, the case D2 > Di. $i(0,i2) > -D if 
and only if I? < D^ . 



Proof. The Jacobian matrix J of ([5]), at point [xi, X2), is given by: 

ds2 



dfi ^dh ^f 
osi 0S2 



D 



Xl 



5/2 ^ 5/2 
OSi OS2 



where the functions are evaluated at (s™ — xi, s™ 
by: 



dh 

dsn 



X2 +f2-D 



Xl — X2) ■ The Jacobian matrix at F^ is given 



J° 







/2(sr,. 



D 



The eigenvalues are Di — D and D2 — D. Thus, if Z? > max(Z?i, D2) then i^° is a stable node. It is 
an unstable node if Z? < min(Z?i, _D2)- It is a saddle point if min(Di, D2) < D < max(Di, Z?2)- O 

The conditions of existence of the boundary equilibria F^ and F^, and their nature, are stated 
in the following lemmas. 

Lemma 2. An equilibrium F^ = (xi,0) exists if and only if D < Di. If it exists then it the unique 
equilibrium on the positive Xi semi-axis. If Di < D2 then F^ is a saddle point for all D < Di. If 
D2 < Di, then F^ is a saddle point for all < D < D^ and a stable node for all D^ < D < Di. 

Proof. An equilibrium F^ = (Si,0) exists if and only if xi — xi e]0, s™[ is a solution of (|6]). Let 
ipiixi) = $i(a;i,0). Then 

^'iixi)^-^isT-xi.sl"+xi) + ^isT-xi,sT + xi). 

OSi OS2 

By assumption H3, V'K^^i) < 0. Since ■01 (0) = Di, and 01(5™) = 0, equation ^ admits a solution 
in the interval ]0, s™[ if and only if 7D < Di. If this condition is satisfied then ^ admits a unique 
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solution since the function ipii-) is decreasing, see Figure[5] The Jacobian matrix at F^ is given by: 



J' = 


r 5/1. a/i. 9/1 _ 1 

Xi + Xl ^ Xl 

osi 0S2 0S2 

f2-D _ 




re evahiated at (s™ — xi,s™ + xi). The eigenvalues are 


Xl,s'^+Xl] 


- D = <^2{xi,0) - D, and - 


T^—Xl + -^—2^1 < 



Thus F^ is a saddle point if $2(a;i, 0) > D. If _Di < _D2, this condition is satisfied for all D < Di. 
If D2 < Di, it is statisfied for al\0 < D < D^, see Figure [H F^ is a stable node ii D3 < D < Di 
andL>2<-Di- □ 

Lemma 3. An equilibrium F^ = (0,a;2) exists if and only if D < D2. If it exists then it the unique 
equilibrium on the positive X2 semi-axis. If D2 < -Di then F^ is a saddle point for all D < Z?2. // 
Di < D2, then F^ is a saddle point for all < D < D4 and a stable node for all D^ < D < D2. 

Proof. An equilibrium F^ — (0,a;2) exists if and only if X2 — X2 e]0, S2"[ is a solution of ([7]). Let 
Mx2) = '^2{0,X2). Then 

i^'2{x2)^-^isT,sT-i2). 



ds2 



By assumption H4, ^'2(2^2) < 0. Since "02(0) = -D2, and ^"2(5™) — 0, equation ([7]) admits a solution 
in the interval ]0, s'!^[ if and only if D < I?2- If this condition is satisfied then (O admits a unique 
solution since the function ?/'2(-) is decreasing, see Figure[3] The Jacobian matrix at F^ is given by: 

fi-D 



J' 



df2. 

OSi 

where the functions are evaluated at {s^{^^s™ 



df2-. 



-X2 



9/2 



-X2 



/i(sr,4" 



ds2 ds2 

~ X2). The eigenvalues are 
9/2- 



X2)- D = $1(0, X2) - D, and - ^X2 < 0. 

OS2 



Thus F^ is a saddle point if $1(0,5:2) > D. If D2 < Di, this condition is satisfied for all D < D2. 
If Di < D2, it is statisfied for all < D < D4, see Figure 121 F"^ is a stable node ii D^ < D < D2 
and Di < D2. □ 



Let us discuss now the conditions of existence of positive equilibria F* , and their number. An 

equilibrium F* — (x*, Xj) exists if and only if Xi = x^, X2 = xj is a solution of ([5]) lying in S. One 

has 

5$i dfi , ,„ „, , 

0x2 OS2 

By assumption H3, this partial derivative is positive. Hence, equation $i(a;i,X2) = D defines a 
function X2 = Fi{xi) such that Fi(a;i) — when D < Di. Recall that xi = xi is the solution of 
Q which, according to Lemma [2] exists and is unique, if and only if Z? < Di. One has 

p,f N _ dxi (Xi,riyxi)) _ Osi + ds2 _ T Osi ^ , 
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Hence the function Fi is increasing. Since $1(5™, 0) = 0, the graph Fi of Fi has no intersection 
with the right boundary of the domain 5, defined by xi — s™. This graph separates S in two 
regions denoted as the left and right sides of Fi, see Figure IH One has also 



a$2 

dx2 



dh 
ds2 



(sT 



Xi,S2 + Xi - X2) 



By assumption H3, this partial derivative is positive. Hence, equation $2(2;!, X2) — D defines a 
function X2 — F2{xi) such that ^2(0) — X2 when Z? < Z?2- Recall that X2 = X2 is the solution of 
([7]) which, according to Lemma |3] exists and is unique, if and only ii D < D2. One has 



F^{x^) 



dxi y-^l^^2yxi)j ^ Q^_^ ^ Ss2 _ 1 dsi 



dx'^ 



■{xi,F2{xi)) 



dh 
ds2 



I 



dh 
dsi 

W2 
ds2 



> 1. 



Hence the function F2 is increasing. Since $2(2^1, s^^+xi) = 0, the graph F2 of F2 has no intersection 
with the top boundary of the domain S, defined by X2 = si^ + xi. Thus the point at the very right 
of F2 lies necessarily on the right boundary of S, defined by xi — s™. Hence it lies on the right 
side of Fi, see Figure H) 






Figure 4. On the left, the left and right sides of Fi. On the center, the point A 
at the very left of F2 lies on left side of Fi: there are generically an odd number of 
intersections (3 in this example). On the right, the point A at the very left of r2 
lies on right side of Fi: there are generically an even number of intersections (2 in 
this example). 



The graphs Fi and r2 can intersect or not, see Figures SI [S] and [SI If they intersect at some 
point F* = {xl,X2) then F* is a positive equilibbrium. If the point A at the very left of F2 lies 
on left side of Fi then Fi and F2 intersect in at least one point F* — {xl,X2). They can have 
multiple intersections. Generically they have an odd number of intersections (see Figure [U center). 
If the point A at the very left of F2 lies on right side of Fi then Fi and F2 can intersect or not. 
Generically they have an even number of intersections (see Figure^ right). The nature of a positive 
equilibrium F* is stated in the following lemmas. 

Lemma 4. If an equilibrium F* = (xl,X2) exists then it is a stable node if F[(xl) > F2{xl). It is 
a saddle point if the opposite inequality is satisfied. 
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Proof. The Jacobian matrix at F* is given by: 



J* 



Jh . 

where the derivatives are evaluated at (s™ 



dfi , 



dh 



dh 






9S2' 
5/2, 

- x*^. Notice that 
dh 



<0 



and 



det(J*) 






F[{A) 



'dhdh_dhdh' 

dsi ds2 ds2 dsi 

By Assumptions H3 and H4, the product of the partial derivatives is negative. Therefore, the 
determinant is positive if F^'(x*) > i^jl^^i) ^^d negative if the opposite inequality is satisfied. 
Hence the equilibrium F* = (x*, 0:2) is a stable node if Fl{xl) > ^2(0;*). It is a saddle point if the 
opposite inequality is satisfied. D 

The number of equilibria of ([5]) and their nature are summarized in the next theorem. 

Theorem 1. (1) If D < nrin(Z)i, 132) then ^ admits the trivial equilibrium F'^ which is an 

unstable node, the boundary equilibria F^ and F^ which are saddle points, and at least one 
positive equilibrium F* . If F* is the unique positive equilibrium then it is a stable node. 
Generically, the system has an odd number of positive equilibria which are alternatively 
stable nodes and saddle points, the one at the very left of these positive equilibria is a stable 
node. 
(2) If iDhi{Di, D2) < D < maii{Di, D2), four subcases must be distinguished 

(a) // Di < D2 and Di < D < D4 then 0j admits the trivial and boundary equilibria F'^ 
and F"^ , which are saddle points and at least one positive equilibrium F* . If F* is the 
unique positive equilibrium then it is a stable node. Generically, the system has an odd 
number of positive equilibria which are alternatively stable nodes and saddle points, the 
one at the very left of these positive equilibria is a stable node. 

(b) If Di < D2 and D4 < D < D2 then ([^ admits the trivial equilibrium F'^ , which is a 
saddle point, and the boundary equilibrium F^ , which is a stable node. Generically, the 
system can have an even number of positive equilibria which are alternatively saddle 
points and stable nodes, the one at the very left, of these positive equilibria is a saddle 
point. 

(c) // D2 < Di and D2 < D < D3 then 0) admits the trivial and boundary equilibria F'^ 
and F^ which are saddle points and at least one positive equilibrium F* . If F* is the 
unique positive equilibrium then it is a stable node. Generically, the system has an odd 
number of positive equilibria which are alternatively stable nodes and saddle points, the 
one at the very left of these positive equilibria is a stable node. 

(d) If D2 < Di and D^ < D < Di then ^B^ admits the trivial equilibrium F^ , which is a 
saddle point, and the boundary equilibrium F^ , which is a stable node. Generically, the 
system can have an even number of positive equilibria which are alternatively saddle 
points and stable nodes, the one at the very left of these positive equilibria is a saddle 
point. 
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(9) 



(3) If D > raaK{Di, D2) then ^ admits the trivial equilibrium F^ which is a stable node. 
Generically, the system can have an even number of positive equilibria which are alterna- 
tively saddle points and stable nodes, the one at the very left of these positive equilibria is 
a saddle point. 

4. Growth functions of Monod type 

In this section we consider growth functions /i and /2 of the following form 

misi . , , m2Si 



/i(si,S2) 



/2(S1,S2) 



{Ki + si)iLi + S2) ^-^ - -' iK2 + si){L2 + S2) 

Such functions are simply the product of a Monod function in si by a decreasing functions of 
S2. Such functions are currently used in biotechnology when the growth of a functional species 
is limited by a substrate while inhibited by another one. Such situations are common in water 
treatment technology like in the denitrification (limited by the nitrate and inhibited by the dissolved 
oxygen) or in the anoxic or anaerobic hydrolysis (limited by the slowly biodegradable substrates 
while inhibited by an excess of oxygen) processes which are modeled this way (cf. [IS])- 

One can readily check that © satisfy Assumptions HI to H4. By straighforward calculations 
one has 

Dxl +[mi + D {Ki -Li+ sf - 4")] - misl" + D {Ki + 4") (Li + s^) 



Fiixi) 



F2{xi) = 



Xl) 



Dxl 



Pl2 



D (L2 - K2 



D{Ki + si" 

4»-sf)] +TO24" 



D (K2 + si 



m2 - D{K2 + s\" - Xl) 



Hence equation Fi(a;i) = F2{xi) giving the abscissa of positive equilibria is an algebraic equation 
of degree 2. Thus, it cannot have more than two solutions. Hence, the situation depicted on the 
center of Figure SI of three positive equilibria, is excluded. However, the situation depicted in the 
right of Figure m with two positive equilibria can occur. For instance, consider the following values 





3/5 < D< 8/9 




F° F' 
8/9 < D< 1 




1< D< 6/5 




6/5 < D 
Figure 5. Relative positions of the isocline ii = (in red) and ±2=0 (in green). 



of the parameters 
(10) mi 

Then 



i, m2 =4, Ki=L2 = 1, Li=K2 = 2, s™ = 4" 



Di = 6/5, D3 = 8/9, D2 = 3/5. 
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There is another bifurcation value, D = 1 which correspond to the case when the graphs Fi and 
r2 are tangent, see Figure [6l For this example five cases can occur, see Figure [5j 



Proposition 2. Consider system {^ where /i and /2 
(1 



are given 



by (^ with parameters I110\). The 



when D < 3/5, the system, has four equilibria, F^ which is an unstable node, F^ and F"^ , 
which are saddle points and F* , which is a stable node. This is case (1) of Theorem]^ with 
a unique positive equilibrium. 

(2) when 3/5 < D < 8/9, the system has three equilibria, F'^ and F^ , which are saddle points 
and F* , which is a stable node. This is case (2.c) of Theorem]^ with a unique positive 
equilibrium. 

(3) when 8/9 < D < 1, the system has four equilibria, i^° and F* , which are saddle points 
and F^ and F2 , which are stable nodes. This is case (2.d) of Theorem]^ with two positive 
equilibria. 

(4) when 1 < D < 6/5, the system has two equilibria, F^ , which is a saddle point and F^ which 
is a stable node. This is case (2.d) of Theorem\^ with no positive equilibrium. 

(5) when D > 6/5, the system has one equilibrium, F^ , which is a stable node. This is case (3) 
of Theorem]^ with no positive equilibrium. 







D = 3/5 D = 8/9 D = 1 D = 6/5 

Figure 6. The non hyperbohc cases. When D — 6/5, F° and F^ coalesce. When 
D = 8/9, Fi and F^ coalesce (saddle node bifurcation). When D — 1, F^ and F2 
coalesce (saddle node bifurcation). When D = 6/5, F° and F^ coalesce. 



In the case when 8/9 < 13 < 1 a bistability phenomenon occurs. According to the initial 
condition, both species can coexist at equilibrium F2 , or species X2 goes to extinction at equilibrium 
F^ . This phenomenon is illustarted numerically with D = 0.95 in Figure [71 

For the following values of the parameters 



(11) 



TOi =8, m2 = 7, Ki = K2 



1, Li = 3/2, 



the bifurcational values are Di — 4/3 and D2 = 21/16. If D > max(Z3i, 1)2), for instance for 
D = 3/2, one obtains a bistability phenomenon corresponding to case (3) of Theorem [U with two 
positive equilibria. According to the initial condition, both species can coexist at equilibrium F2, 
or both species go to extinction at equilibrium F'^. This phenomenon is illustarted numerically in 
Figure [H 
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F"F^ 



Figure 7. Numerical solutions in the bistability case D — 0.95 and parameters 
values (flOl) . On the left, the separatrix (in green) of the saddle point F^ separate 
the domain S in two region which are the basins of attraction of the boundary 
equilibrium point F^ and the positive equilibrium point i^j*- On the center, the 
phase portrait. On the right, the isoclines. 






Figure 8. Numerical solutions in the bistability case D — 1.5 and parameters 
values pT|) . On the left, the separatrix (in green) of the saddle point _F\* separate 
the domain S in two region which are the basins of attraction of the boundary 
equilibrium point F'^ and the positive equilibrium point F2 . On the center, the 
phase portrait. On the right, the isoclines. 



5. Global analysis 

Let us establish first that ([5]) admits no periodic orbit nor polycycle inside S 

Theorem 2. There are no periodic orbits nor poly cycles inside S. 

Proof. . Consider a trajectory of ([5]) belonging to S. Let us transform the system ([S|) through the 
change of variables ^1 = ln(a;i), ^2 = ln(x2). Then one obtains the following system : 



(12) 
We have 



6 = /ii (a , 6 ) := /2 (4" - e«^ , 4" + e«^ - e«^ ) - D. 
dhi dh2 e^ dfi , dfi . a/2 

0^1 0?2 OSi dS2 OS2 
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From Dulac criterion [IT], we deduce that the system ([T^ has no periodic trajectory. Hence ((S]) 
has no periodic orbit in S. D 

Theorem 3. Assume that system 0j has at m,ost one positive quilibrium F* , then for every initial 
condition in S, the trajectories of system ^ converge asymptotically to : 

• F* if D <Tnva{Di,D2). 

• F* if Di < D2 and Di < D < D^ 

• F'^ if Di < D2 and D^ < D < D2. 

• F* if D2 < Di and D2 < D < D3. 

• Fi ifD2 < Di and D3 < D < Di. 

• F° ifui'Ayi{Di,Di) < D. 

Proof. We restrict the proof to the situation where D < inui{Di,D2)- The other cases can be done 
similarly. Let2;i(0) > 0,a;2(0) > and w the w-limit set of (a;i(0), 0:2(0)). w is an invariant compact 
set and to C S. Assume that w contains a point M on the xia;2 axis : 

• M can't be F° because F° is an unstable node and can't be a part of the w-limit set of 
(xi(0), 2:2(0)), 

• If M e]a;i,s™] X {0} (respectively M e {0}x]a;2, Ss"])- As uj is invariant then -f{M) C uj 
which is impossible because w is bounded and 7(Af) —]xi, +oo[x{0} (respectively 7(A/) = 
{0}x]x2,+oo[), 

• If M g]0,xi[x{0} (respectively M e {0} , x]0,X2[). w contains 7 (M) =]0,xi[x{0} (respec- 
tively j{M) = {0}x]0,X2[)- As oj is a compact, then it contains the adherence of 7(M), 
[0,xi] x {0} (respectively {0} x [0,iE2])- In particular, w contains F'^ which is impossible, 

• If M = i^^ (respectively M = F"^). w is not reduced to F^ (respectively to F^). By Butler- 
McGehee theorem, to contains a point P of (0, +00) x {0} other that F^ (respectively of 
{0} X (0, +00) other that F^) which is impossible. 

Finally, the w-limit set don't contain any point on the xia;2 axis. System ([5]) has no periodic 
orbit inside S. Using the Poincar-Bendixon Theorem [TTi, F* is a globally asymptotically stable 
equilibrium point for system ([5]). D 

Theorem 4. Assume that system ^Bj) has at most one positive quilibrium F* , then for every initial 
condition in M.\_ , the trajectories of system (0j converge asymptotically to: 

• E* if D < min{Di,Di). 

• E* if Di < D2 and Di < D < D4. 

• E"^ if Di < D2 and D4, < D < D2. 

• E* if D2 < Di and D2 < D < D3. 

• E'^ if D2 < Di and D3 < D < Di. 

• E° i/max(Di,D2) < D. 

Proof. Let {si{t) , xi{t) , S2{t) , X2{t)) be a solution of 12]). From ([3]) and (H]) we deduce that 

si(i) = sr-a;i(t) + ifie-^* and S2(i) = 4" + a:i(t) - X2(i) + ^^26"^*, 

where Ki = si(0) + xi(0) - sf and K2 = S2(0) + 2:2(0) - a;i(0) - s|". Hence {xi{t),X2{t)) is a 
solution of the nonautonomous system of two differential equations : 

ii - [fi{sT-x,+K,e-''\s'J'+x,-X2 + K2e-^')-D]xi, 
(13) 

±2 = [f2{sT-xi+Kie-^',s'2"+^i-^2 + K2e-''')~D]x2. 



This is an asymptotically autonomous differential system which converge to the autonomous system 
(O. The set fl is attractor of all trajectories in R^ and the phase portrait of system reduced to 
il (O contains only locally stable nodes, unstable nodes, saddle points and no trajectory joining 
two saddle points. Thus we can apply Thiemes's results [12] and conclude that the asymptotic 
behaviour of the solution of the complete system ([T3|) is the same that the asymptotic behaviour 
described for the reduced system ^ and the main result is then deduced. D 

6. The anaerobic digestion process : An example of a synthrophic relationship 

"Methane fermentation" or "anaerobic digestion" is a process that converts organic matter into 
a gaseous mixture mainly composed of methane and carbon dioxide (CH4 and CO2) through the 
action of a complex bacterial ecosystem (cf. FiglH]). It is often used for the treatment of concen- 
trated wastewaters or to stabilize the excess sludge produced in wastewater treatment plants into 
more stable products. There is also considerable interest in plant-biomass-fed digesters, since the 
produced methane can be valorized as a source of energy. It is usually considered that a number of 
metabolic groups of bacteria are involved sequencially. 

One specific characteristic of the anaerobic process is that within such groups, there exists pop- 
ulations exhibiting obligatory mutualistic relationships. Such a syntrophic relationship is necessary 
for the biological reactions to be thermodynamically possible. In the first steps of the reactions 
(called "acidogenesis"), some hydrogen is produced. In El Hajji et al.[5], this production of hydro- 
gen at this reaction step was neglected (compare Fig|9]with Fig.l of [5]). This hypothesis constitue 
the first novelty with respect to [5^ . It is to be noticed that an excess of hydrogen in the medium 
inhibits the growth of another bacterial group called "acetogenic bacteria" . Their association with 
H2 consuming bacteria is thus necessary for the second step of the reaction to be fulfilled. Such a 
syntrophic relationship has been pointed out in a number of experimental works (cf. for instance 
the seminal work by [3]). Let us consider the subsystem of the anaerobic system where the VFA (for 
Volatile Fatty Acids) are transformed into H2, CH4 and C02. We can formalize the corresponding 
biological reactions as a first bacterial consortium Xi (the acetogens) transforming ^i (the VFA) 
into 5*2 (the hydrogen) and acetate (cf. Fig[S]). Then, a second species X2 (the hydrogenotrophic- 
methanogenic bacteria) grows on 82- In practice, acetogens are inhibited by an excess of hydrogen 
and methanogens by an excess of VFA. Thus, it is further assumed that Xi is inhibitied by S2 and 
X2 by Si. The last inhibition relationship constiute the second novelty with respect to f5^. This 
situation is precisely the one considered within the model p]). 

We have proposed a mathematical model involving a syntrophic relationship of two bacteria. 
It results from this analysis that, under general and natural assumptions of monotonicity on the 
functional responses, the stable asymptotic coexistence of the two bacteria is possible. 
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